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Abstract 

The issue of different parameterizations of the axisymmetric vesicle 
shape addressed by Hu Jian-Guo and Ou-Yang Zhong-Can [ Phys.Rev. 
E 47 (1993) 461 ] is reassesed, especiaUy as it transpires through the 
corresponding Euler - Lagrange equations of the associated elastic en- 
ergy functional. It is argued that for regular, smooth contours of vesi- 
cles with spherical topology, different parameterizations of the surface 
are equivalent and that the corresponding Euler - Lagrange equations 
are in essence the same. If, however, one allows for discontinuous 
(higher) derivatives of the contour line at the pole, the differently 
parameterized Euler - Lagrange equations cease to be equivalent and 
describe different physical problems. It nevertheless appears to be true 
that the elastic energy corresponding to smooth contours remains a 
global minimum. 

*0n leave from: Department of Theoretical Physics, J.Stefan Institute, P.O.Box 100, 
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Recently Hu Jian-Guo and Ou-Yang Zhong-Can [Q] argued to the effect 
that different parameterizations of the axisymmetric shape of a hpid vesi- 
cle lead to different equations governing its shape. Should this assertion be 
regarded as generally valid, it leads to disturbing ambiguities in the results 
of the computation of equilibrium shapes of the elastic vesicles. Also many 
results derived till now [0, |^ and based on the assumption that the parame- 
terization of the surface does not alter its equilibrium form, would be cast in 
grave doubt. It is thus of utmost importance to establish the precise limits of 
the claims made in Ref.l, particularly since it is quite clear from the general 
theory of the variational calculus 0] that these claims can not be generally 
valid. 

A criticism of the variational methods employed in [|I| has already been 
voiced by Jiilicher and Seifert 0, that based their arguments on the incor- 
rect way the boundary terms were treated in the derivation of the variational 
equations for general topologies. Limiting ourselves to the spherical topol- 
ogy we will, in what follows, investigate the variational shape equations for 
different parameterizations and their range of validity. 

Though all our formal derivations and associated discussion will be cen- 
tered on the bilayer - couple model, the modifications brought to the formal- 
ism in the framework of the spontaneous curvature model - the two are just 
different limiting cases of the generalized bilayer - couple model - are known 
[^] and do not in any respect modify the conclusions reached below. 

The formulation of the problem treated here is simple. We start with the 
variational principle for an elastic vesicle under the constraint of constant 
volume, V, constant surface area. A, and constant difference between the 
areas of the two membrane surfaces in contact, A A. The equilibrium shape 
of the vesicle is obtained from the first variation 



where the bending energy can be expanded in a series with respect to the 
principal curvatures ci and C2 of which the first term is of the form 



In the above equation we omitted the Gaussian elastic term, as we consider 
only vesicles of fixed (spherical) topology. Furthermore the elastic energy, 
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area, volume and area difference are all normalized with respect to the cor- 
responding values for the spherical shape 0], while L, M and are the 
Lagrange multipliers corresponding to the constraints of fixed area, volume 
and area difference. Furthermore da is the normalized area element of an 
axisymmetric vesicle, i.e. j da = 1. 

We start by choosing the relevant variables describing the shape of an 
axisymmetric vesicle. These are the angle il){x) between the tangent to the 
contour and the x - axis, and the contour line itself z{x). The z - axis 
coincides with the rotational axis and is perpendicular to the The 
connection between the contour line and il){x) is 

= -tanV'(x), (3) 

dx 

while the two principal curvatures of an axisymmetric shape have the form 

sin?/'(a;) dip{x) , , 

Cp[x) = and Cm{x) = cosip{x) — - — . (4j 

The volume and the surface area difference can be obtained in the form 

V = ^ <j) x^ Cp da 
AA = I ^ (cp + Cm) da. (5) 

We now analyse four different derivations of the variational shape equa- 
tions , all stemming from Eq.|l], of which three are based on axial radius, area 
and contour arclength parameterizations, while the last one is a "general" 
shape equation 0, that for an axisymmetric shape reduces to j- parame- 
terization. To the analysis of ^ that contains the comparison of the first 
and the third of these parameterizations and the "general" shape equation 
we thus add an additional derivation of the shape equation, where the area 
itself is taken ClS db variational parameter. 

If we start by taking the axial radius x as the parameter of the shape, we 
obtain the contour line z{x) from the minimization of the action functional 
C which can be cast into a dimensionless form 

C = lJLJx,^P{x),^)dx, (6) 
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The Euler - Lagrange equation corresponding to the x - parameterization can 
be obtained from Eq.^ in the form: 

2 
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One notes here that this equation is of the second order with respect to 
the derivatives of ip{x). The contour hne follows after a solution of Eq.^ is 
inserted into Eq.^ 

We now proceed with the surface area parameterization of the vesicle 
shape. We take the axisymmetric area element in the form 



da 



x doc 
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together with %p = ip{a). In this parameterization the Lagrangian is obtained 

dip (a) 



as 



da 
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with 



X, ip{a) 



dip [a) 
da 



f X dip sin ip\ ^ , ^ . , X dip sinip\ 

= — r- + -] -L-Mxsmip-Nl — J- + + 

\2 da X J \2 da x J 



The last term in the above equation is due to the constraint Eq.|^ in the 
form of a differential relation between the surface area element and the axial 
radius. 
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In the a - parameterization we now obtain two Euler - Lagrange equations, 
one for ip field, which is of the second order in the derivatives of ip, and the 
other one for 7q, field, which is of the first order. Combining these two 
equations in analogous way as in Ref.l into a single equivalent equation for 
the ip field of a higher (third) order , we are led to 

^ZM = o ^ n{x)= const., (12) 
ax 

which is now the Euler - Lagrange equation in the a - parameterization. The 
7i(a;) function is the same as has been defined in Eq.|^. 

Taking the arclength of the vesicle contour as a parameter, thus 

ds = — ^JT^, (13) 
and = ipls), we obtain the Lagrangian in the form 

C=ULJx,^Pis),^]ds, (14) 



ds 



where 



lJxMs),^] = x(^ + '^X -Lx-Mx^sin^-Nx(^ + '^] + 
\ ds J \ds X J \ds X J 

+ 7s - cos?/'^ . (15) 

The last term in the above equation again ensures that the geometric con- 
straint Eq.|13| is obeyed along the whole of the contour. 

The Euler - Lagrange equations are once again two, one for the field 
and the other one for the 7^ field. Combining the two of them into a single 
higher order equation for ip we derive 

^ + ^fl-^Uo - ^ = const., (16) 
dx X \ Cp{x) J Cp{x) 

where again Cp and Cm are the two principal curvatures and the function 7i(a;) 
has again been defined in Eq.& 
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The next alternative derivation of the Euler - Lagrange equations for 
the shape of the vesicle proceeds by describing the variation of the contour 
line in terms of infinitesimal deformations of the shape in the direction of 
the local normal to the unperturbed shape [0. For an axisymmetric shape 
this is equivalent to taking — as a parameter describing the shape. This 
procedure leads yet to a new Euler - Lagrange equation, that in the case of 
an axisymmetric vesicle reduces to Eq.(7) of Ref.l. After some algebra we 
were able to write the corresponding Euler - Lagrange equation in terms of 
7^(x), that in this case assumes the form 

dHix) H(x) „ X /..„s 

— ^ + = n(x) X = const. (17 

ax X 

The same equation has been recently derived by a somewhat different pro- 
cedure also by Zheng and Liu [§] in their analysis of the axisymmetric shape 
equations. 

We were thus able to show that for these four different derivations of the 
Euler - Lagrange equations, Eqs.^|l2],|l6| and giving the contour line of 



the axisymmetric shape, invariably reduce to a general equation of the form 
7Y(x) f{x) = const., where the form of the function f{x) and the value of 
the constant depend on the type of derivation. 

Let us now investigate how the solutions of the different Euler - Lagrange 
equations differ between themselves. This is most easily seen if we investigate 
the nature of the contour line close to the pole, i.e. x— >^0,s— >-0, a— >-0 
in different parameterizations. By expanding the different Euler - Lagrange 
equations in the vicinity of the pole we get the following limiting behavior. 
In the X - parameterization the expansion of Eq.H leads to 

Cp{x) =a + bx^ + (18) 

where the two constants are determined from the constancy of volume (V), 
surface area (A) and the surface area difference {AA in the bilayer - couple 
model). In the area parameterization the corresponding limiting behavior 
following from Eq.|l^ is 

Cp{x) = a + b\x\ + . . . , (19) 

where the constant a is determined from the values of V, A, AA, while b = 
const./ 3 where the constant is the same as in Eq.O. The arclength param- 



6 



eterization leads to the following limiting law 



a + b\x\ + ..., (20) 



with — = const., where the constant is the same as in Eq. ITB 



Finally for the Euler - Lagrange equation Eq.|T^, corresponding to locally 
normal shape variations introduced by Ou-Yang and Helfrich [0 , one gets 

Cp{x) = a log |x| + . . . , (21) 

with a = const./ 2, where the constant is the same as in Eq.|l^. Interestingly 
enough an equation of the form Eq. |2T| has been derived recently by Naito 
et al. as a general solution (not just a limiting form as is the case in Eq. 



21|) of their shape equation, which is just the Euler - Lagrange equation in 
the Ou-Yang - Helfrich parameterization. 

One can easily demonstrate that the lowest power n in the expansion 
Cp (x) = a + 6x" + . . . compatible with at least one of the four Euler - Lagrange 
equations is = "0" (in the sense that "x"^'^" = logx ), corresponding to 
the Ou-Yang - Helfrich parameterization. On the other hand if we choose 
the constants in Eqs.|T2|, |16| and |17| as zero, all the Cp's reduce to the limiting 



form Eq. 

The contour line close to the poles is then obtained in an approximate 
form as 

'^^^ = - tan ^(x) = - (xcp(x)) . (22) 
ax 

It is thus clear that only the form Eq.^ keeps the contour line and all its 
higher derivatives finite close to the pole. The other limiting forms lead to 
discontinuities in the third derivative of the form 

^ ^ -2b sgn(x) (23) 

for the limiting laws Eqs.^ and |2D|, or of the form 

d^z{x) ^ a 



dx^ 



\x\ 



sgn{x) (24) 



for the limiting case Eq.pl|. 

The existence of discontinuities in (higher) derivatives of the contour z{x) 
is in general connected with point source terms in the Lagrangian. Let us 
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investigate this connection in the case of the Ou-Yang - Helfrich parameter- 
ization. An external force would in general introduce a term of the form 



into the energy of the vesicle, where f{x) is the linear density of the force 
in the direction of z - axis and the numeric coefficient has been chosen for 
later convenience (if there were external momenta acting on the contour we 
would in general have to include those contributions too). As Eq. refers 
to a discontinuity at the pole, we presume that the force density has the form 
f{x) = /o5(x), where /o is a constant. Deriving now the Euler - Lagrange 
equation for the ansatz Eq.^ we can cast it in the familiar form 



force acting at both poles. Instead of using the language of point sources in 
the Lagrangian we could also invoke external constraints: if the separation 
between the poles is kept constant Az = const., then /o is just the conjugate 
variable to Az. 

Either way the conclusion is that the non-analyticities of the Hu - Ou- 
Yang type in the derivatives of the contour line at the poles of the form Eq.p^ 
exist only if there are external forces acting on the poles, pulling (pushing) 
them apart (together). More complicated external sources (constraints) are 
needed in the case of other parameterizations. 

Another conclusion pertinent to the above discussion is that no matter 
what the magnitude of /o might be, the corresponding membrane bending 
energy will always be larger than in the /o = case. On Fig.l we present 
a result of numerical computation of J-'b Eq.^ for a range of /o values while 
keeping V and AA constant ( V = 0.95 and AA = 1.0129). Obviously /o = 
case represents a global minimum of this energy. 

One is thus lead to the following conclusion regarding the different forms 
of the Euler - Lagrange equation, stemming from the various parameteriza- 
tions of the vesicle shape. If one demands that the contour line be a smooth 
function, without any discontinuities in its values or its derivatives, then all 
the different Euler - Lagrange equations reduce to a single equation of the 
form T-C{x) = 0, identical to the x - parameterization result. This is con- 
sistent with the parameterization invariance of the variational problem, but 




(25) 
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is at odds with the generahty of the claims made in Ref.l. On the other 
hand if one allows for discontinuities in the (higher) derivatives of the con- 
tour line of the vesicle, then each Euler - Lagrange equation corresponding 
to different parameterizations gives different equilibrium shapes. However, 
these shapes correspond to different physical problems. It appears that only 
the parameterization where 7i(x) = describes a completely unconstrained 
vesicle of spherical topology. The other parameterizations describe the shape 
of a vesicle under external forces and are not relevant to the problem of the 
shapes of an unconstrained vesicle, but might provide additional insight into 
the shapes of constrained vesicles. 

The conclusion of this paper, contrary to the assesments of is that 
for "smooth" axisymmetric shapes the way in which one derives the shape 
equations matters little. All the different parameterizations employed in the 
literature are equally correct and their corresponding variational equations 
are parameterization invariant - a results of the general theory of variational 
calculus of long standing |Q. Should one, however, want to include 'patho- 
logical' shapes with various discontinuities in the higher derivatives of the 
contour line one should first of all start with a more elaborate form of the 
curvature energy including the external forces or geometrical constraints that 
lead in each particular case to these discontinuities. For these cases no general 
form of the free energy is, however, presently in sight and the discussion of 
the variational equations for those cases should reflect the explicit constraints 
for each particular case. 
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Figure Captions 

Fig.l: Membrane bending energy Th obtained from numerical integration of 
Eq.^ for different values of the force /q. The normalized volume is ^ = 0.95 
and the normalized area difference is = 1.0129. 
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